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I. INTRODUCTION 


In the quest of an UV-complete quantum field theory of gravity, the search for a fixed point 
using functional renormalization group methods has reached the point where one may hope to 
go beyond finitely many couplings and study entire functional classes of truncations. The best 
studied case is that of f{R) actions, where a fixed point is known to exist, and to exhibit nice 
stability properties, when / is a polynomial m- The most advanced calculations have now 
reached order ISIS]. However, the radius of convergence of the Taylor series of / around 
the origin is finite and there is not much to be gained by pushing the expansion much further. 
Rather, one would like to find a scaling solution for the whole function /. Several studies 
have shed light on various aspects of this issue but have so far failed to reach a convincing 
conclusion, at least in four dimensions iHII]. An important fact that has been pointed out 
in [8] is that the equation of mm does not admit complete solutions. The simpler equation 
proposed in |6] admits solutions at least for positive R but then it was shown in m that all 
perturbations around them are redundant, i.e. can be absorbed by field redefinitions. One thus 
has to find a “better” equation, i.e. one admitting a discrete set of solutions with non-redundant 
perturbations, or else show that no such equation exists. In order to gain some understanding 
of what may be wrong with the equations of laElE], it has been shown in [T^] that the use of 
background-dependent regulators in the flow equation for a scalar field can artificially lead to 
similar pathologies. It is therefore important to understand whether different ways of applying 
the background field method could solve this issue. 

In this paper we will discuss similar problems but in a different context, namely a scalar field 
non-minimally coupled to gravity. We will consider Effective Average Actions (coarse-grained 
effective actions depending on a cutoff k, usually abridged EAA) of the functional form: 

rfc[<^, g] = j d<^x ^ (vi^) - Fi<p)R + + Sgf + Sgh , (LI) 

where Sgf and Sgh are gauge-fixing and ghost terms. The usual Einstein-Hilbert action is 
contained in this truncation as the constant (^-independent) part of the action, while switching 
off gravity (i.e. setting g^u = 6fj,u) reduces the system to the well-studied Local Potential 
Approximation (LPA) of the scalar field. 

There are several good reasons to study such actions. On one hand they may have direct 
applications to cosmology m- At the classical level, they are related via some field redefinitions 
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to the f{R) class of actions. ^ Whether the classical equivalence can be preserved at the quantum 
level is doubtful at present, but not completely settled. From the point of view of this work, 
they have the important advantage that the scalar and gravity subsectors in isolation are well- 
understood. In particular, in d = 3 the LPA admits a solutions that is a good approximation 
to the Wilson-Fisher fixed point da IS] and pure gravity admits a hxed point that is relatively 
weakly coupled, independent of the presence or absence of higher derivative terms [21H25| . These 
two facts lead us to suspect that the coupled system also should have a fixed point. 


An early study of one-loop divergences of the form (I.l), in fact including also a prefactor 
Z{(f>) for the kinetic term, was made in dSl) see also [20]. Here we shall derive the RG flow of 
14 and Fk from the functional RG equation daim, 


hm = ^STr 




( 1 . 2 ) 


where the dot stands for the partial derivative with respect to t = logk/ko, and Rk is the 
operator which realises the coarse-graining procedure. The flow of I 4 and i4 is a system of 
coupled PDF’s while the fixed-point equations form a system of coupled ODE’s. Solving such 
equations will be the main challenge of this work. 

The flow equations for the theory 0 had been derived earlier in |26j . Fixed point solutions 
had been found for polynomial truncations but they do not have the desired properties, as we 
shall recall in some detail below. It is possible that, just as in the f{R) case, this is due to 
the way the exact equation has been approximated, and in particular to the implementation 
of the background field method. Therefore we derive alternative equations based on a different 
definition of the quantum-background split and a different gauge choice. We will then see that 
the new equations admit nontrivial scaling solutions. 

This paper is organized as follows. In section II we discuss the issues encountered by the old 
flow equations. In section III we motivate the use of the exponential parametrization for the 
metric and of “physical” gauge choices, in particular of the “unimodular” gauge. We then derive 
new flow equations for F and V, valid in any dimension. For simplicity we discuss hrst the case 
when the terms proportional to F in the r.h.s. of the equations are neglected. In sections IV 
and V we discuss in some detail the solutions of these equations in d = 3 and 4 respectively. 
In section VI we briefly describe some results when the terms proportional to F are retained. 


^ This has motivated the study of the RG flow of Brans-Dicke-type actions in m- 
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Section VII contains our conclusions. Three appendices contain a discussion of a multiplicative 
background field method, of functional Jacobians and some results in arbitrary dimensions. 


II. OLD EQUATIONS AND THEIR AILMENTS. 


Flow equations for the functions F and V have been derived in |26) . see also m- They were 
then further simplified by Taylor expanding F and V around (/> = 0 to some finite order, and fixed 
points have been searched within the resulting finite dimensional theory space. In d = 4 the only 
nontrivial solution had constant / and v. It represents a non-interacting scalar field minimally 
coupled to the well-known fixed point of pure gravity in the Einstein-Hilbert truncation. The 
absence of other solutions was perhaps not too surprising, given that such solutions do not 
exist for the pure scalar theory. In d = 3, however, pure scalar theory admits a nontrivial 
scaling solution, the well-known Wilson-Fisher fixed point. In the simplest approximation. 


known as the Local Potential Approximation (LPA), (1.2) reduces to the following equation for 
the dimensionless potential v{(p): 


V = —3v + 



1 

67r^(l -|- v") 


(II.l) 


The solution to this equation can be obtained by a variety of semi-analytic and numerical 
methods. In view of this, there is perhaps greater reason to expect that a nontrivial scaling 
solution may exist also for the system of the scalar field coupled to gravity. However, in [26] 
no such solution was found. A Taylor expansion around (j) = 0 yielded fixed points all of whose 
Taylor coefficients are negative. Even if this corresponded to a genuine fixed point, it would 
not be physically acceptable. If a fixed point existed and was analytic at (/> = 0 it would show 
up as a fixed point for the Taylor coefficients, so the failure to find a fixed point for the latter 
implies that the functional equations do not have a global solution either. The question then 
arises whether this reflects a genuine physical property of the system, or some problem with the 
equations. 

In order to discuss this we will not need to consider the whole equations, it will suffice to 
look at one term that comes from the contribution of the spin two excitations, namely 


V = 


1 

Svr^ 


/ 


f-V 


+ .. 


(II.2) 


In a polynomial expansion of the solution around (/? = 0, it turns out that u(0) < /(O) [26] . 
On the other hand, for large if one expects the solutions to behave like v = and 
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f = Bip'^ + ..., where the dots stand for inverse powers of A solution with these boundary 
conditions would have to cross the singularity at v = f. Although this cannot be ruled out, it 
is likely that the failure to find physically acceptable polynomial fixed points is related to the 
existence of this singularity. 

This conclusion is reinforced by the following two observations. First, that this issue concerns 
the behavior of the dimensionless potential when the dimensionless field = 4> becomes 

large. For d > 2 this is therefore an infrared issue. Second, when F and V are constant, one 
can identify 


V = 


2A 

IGttG ’ 


F = 


1 

IGttG 


(n.3) 


and the fraction //(/ —u) reduces to 1/(1 —2A/k^). The singularity we are discussing is there¬ 
fore a generalization of the well-known infrared singularity at A = k^/2 that appears in most 
treatments of the gravitational flow equation. 

This singularity is an artifact of the way the beta functions for v and / are constructed. The 
inverse propagator for the transverse, traceless spin-2 components is given by: 


F -V^ + 


2 , — 3d + 4 


d(d-l) 


R] -V. 


(n.4) 


The last term comes from the expansion of the in the potential term, which in the standard 
linear background field expansion contains terms of the form The origin of the trouble¬ 

some term in the flow equation for the potential is this propagator, with R set equal to zero. 
But if y / 0, putting R to zero means that we evaluate the flow equation on a configuration 
that is far off shell. Let us see what would happen if we evaluated the equation on shell. For 
constant ([) the trace of the equation of motion of the metric implies 


FR = 


d 

d-2 


V . 


(II.5) 


If we use this relation to eliminate the R term, the spin-2 inverse propagator becomes 

^ ^ (d-l)id-2) 

This would contribute to the flow equation of u a term 


V = 


1 


/ 


/ + 7Fr^ 


+ ... 


(II.6) 


(11.7) 


(d-l){d-2) 

where the troublesome singularity at u = / is no longer present. This is strong evidence that 
the singularity at u = / is unphysical. 
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From the discussion above it is tempting to try and expand the flow equation around a 
solution that is (nearly) on shell. The virtues of such an approach have been discussed previously 
by Benedetti m and Falls [28]. In the following we will not pursue this idea, but rather we 
will employ a different parametrization of the field and choice of gauge fixing that automatically 
avoid the issue. 


III. THE NEW FLOW EQUATIONS 
A. Exponential parametrization 

Instead of the traditional linear quantum-background split g^i, = + h^y we shall use in 

this paper an exponential parametrization 

<7/.. = -g^pie^y (III.l) 

where 5 is a fixed but arbitrary background. This expansion has been used previously in [29|. See 
also m for a recent discussion in a context that is closer to the present one. Some geometrical 
motivation for the use of this formula is given in Appendix A. We assume in this paper that the 
path integral measure is simple when expressed in terms of the field h thus defined. We discuss 
in appendix B the Jacobian relating this measure to the one of the linear parametrization. 

We will use the background metric g to raise and lower indices. Then due to the symmetry 
of gfj^y and gfj^y also the tensor h^y = g^ph^y is symmetric. We have 

g^U — gpy T h^y -\- —h^xh y . (III.2) 

^ ^ (ni.3) 

In contrast to the usual linear split, here also the covariant metric is nonpolynomial in the 
quantum field h^y. Another significant difference is that, due to the formula dete^ = only 
the trace part of h enters in the definition of the determinant, at all orders. As a result 
does not contribute to the action of traceless fluctuations, which are therefore independent of 
the potential. We can split 

(III.4) 

where tr/i = 2dio and is tracefree. Then 

y/g = -|- dw -|- + .. .^ . (III.5) 
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For the scalar field we also expand around a background (p: 

(j) = (p + 54) . 


(IIL6) 


We then expand the action (I.l) to second order in h and 54- Collecting all the terms we hnd 

- ^(trh)V^V,h'^" + ^(tr/i)V2(tr/i) 


-\R^.puah^^''hP^ + - \R{iThf) 


-F'{4) - V2(trh) - + ]^R (tr/i) ) 54 

+ \54{-V^ + V"{4) - F"{4)R)H + \v'{4){t^h)54 + ^^(^)(trh)2 


(IIL7) 


This is identical to equation (6) in [26], which was derived using a linear split, except for two 
terms that are missing here: 


- ]^F{4)R^^^h^phPu - |(^('^) - F{4)R)hp,hP^ . 


(IIL8) 


The latter came from the expansion to second order of the square root of the determinant of g. 
It is absent here because in the exponential parametrization the determinant depends only on 
the trace part of h. 

We then proceed with the York decomposition for the tracefree part of h: 


pv = , 


flU 


(III.9) 


where = 0 and = 0. As usual it is convenient to further redefine 


C = \l■. 


a 


= \/^W-V2 - 


R 


d-l 


a . 


(III.IO) 


Collecting all terms we can rewrite the quadratic action in terms of the independent fields 

fpPT^ ^ ^ 


dxy/g 


1 


Tt 




-V^T 


2R 


d{d-l) 


h 


’j'j'pLu {d — 1)(^ — 2) ! 


Ad? 


a 


(-V2) 


2^ a' 


(-V2) (^-V2 - a'-{d- l){d - 2)uj (^-V2 + I a; 


2{d-l) 


-F'{4)^S4l\l{ 


_V2) ( -V2 - V' + 2d f-V2 + w 




2{d-l) 


+ ^<5<^(_V2 + v"{4) - F"{4)R)54 + V'{4)duj54 + 


(III.ll) 
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Note the absence of from the expansion. Also note that the kinetic operator of the u field is 
not the conformal scalar operator (which has a factor 4 instead of 2 in the denominator). 


B. Gauge choice 


At this point we have to choose a gauge. In order to simplify the equations as much as possible 
we will choose a “physical” gauge, which amounts to putting the gauge-variant components of 
to zero. Such gauges have been discussed earlier in a similar context in [23], see also [3TII33] . 
The transformation of the metric under an infinitesimal diffeomorphism e is given by the Lie 
derivative 


^e9fj.u — T 9fj.p9u^^ T 9updp^^ ■ 


(III.12) 


As usual, we have to define transformations of g and h that, used in (III.l), yield (III.12). The 


simplest one is the background transformation. If we treat g and h as tensors under 5^, i.e. 

= C,gp, ; . (III.13) 

then also 

5[^\e^Y, = C,{e^Y, (111.14) 


and (III. 12) follows. By dehnition, the “quantum” gauge transformation of h is such as to 


reproduce (III.12) when g is held hxed: 

^Pgpu = 0 ; gpp6P {e^Yu = C.g^, . (III.15) 

From the properties of the Lie derivative we have 

C,gp, = C,gpp{Y^Y, + gppCYe^Yu = (V^e^ + V^ep)(e'^)^ + gf^xie-^^C^ie’^Yu ( 111 . 16 ) 
Then we find 


{e-^dPe’^Yu = (e-^CPYu + {epYpPe. + V.eP){eY% (IIL17) 

Expanding for small h we hnd: 

5Ph% = {C,gY. + + [C.g, hY. + 0{ehY ■ (IIL18) 

We note that the first two terms coincide with the quantum transformation when one uses the 
linear background decomposition. In the following we shall only be interested in the functional 












Tk{h]g) for /i = 0. It is therefore sufficient to consider only the hrst term in (III.18). Since 


g = 0 we can write 


V = + 0{h) 


(III.19) 


Using the background g we can decompose the transformation parameter in its longitudinal 
and transverse parts: 


ei^ = + V'"- 


1 


rV'; 


= 0 . 


(III.20) 


The inverse square root of the background Laplacian has been inserted conventionally in the 
dehnition of ip so that it has the same dimension as e^. We can then calculate the separate 
transformation properties of the York-decomposed metric under longitudinal and transverse 
infinitesimal diffeomorphisms. We have 

1 . 2 


M — fTfi . 


= e 




(III.21) 


all other transformations being zero. Note that a and oj are gauge-variant but the combination 


2uj — gV^iT is invariant. In terms of the redefined variables (III.10) we have 




<5^a' = 2W-V2- 


R 


d-l 


ip . 


(III.22) 


First we pick the “unimodular gauge” det^ = det g. ^ In unimodular gravity this is imposed 
as an a priori condition on the metric: by definition the path integral is then over metrics 
with fixed determinant. Here we start from the usual path integral over all metrics and take 
det g = det g as a partial gauge condition. This means that we have to take into account a ghost 
term. To hnd the ghost operator we hrst observe that in the exponential parametrization the 
unimodular gauge condition is 


w = 0 . 


(III.23) 


From (III.21) one then hnds that the path integral must contain a ghost determinant 
det(\/—V^) = det(—V^). As usual, this can be rewritten as a path integral over a real 
anticommuting scalar ghost 

V = / d^xVSc{-V^)c (III.24) 


^ One often just sets det g = 1. This is incompatible with the choice of a dimensionful metric, that we prefer. Also 
note that on a compact manifold the gauge group does not allow one to make constant rescalings of the metric, 
so that the overall scale of the metric remains a physical degree of freedom. We ignore it in the following, since 
it does not affect the running of the terms in our truncation. 
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(This gauge condition has been discussed previously in jMj. There, the ghost was a complex 


scalar. This difference is due to the different definition of if: in (III.20).) 


The unimodular gauge condition completely breaks the invariance under longitudinal in¬ 
finitesimal diffeomorphisms, but leaves a residual gauge freedom that consists of the volume¬ 


preserving diffeomorphisms, which are generated by the transverse vector e^. From (III.22) we 
see that this residual freedom can be fixed by further choosing 


C = o, 


(III.25) 


which gives rise to a ghost determinant det ( \/—“ tI ) = •\/det (— ~ Again, this 


can be written as a path integral over an anticommuting real transverse vector 




(III.26) 


Equations ( III.23|III.25 ) define the “unimodular physical gauge”, which is the gauge condition 
that, unless otherwise stated, will be used in the rest of the paper. 

Before proceeding it is instructive, however, to think for a moment of an alternative choice. 
Since the combination 2oj — is gauge invariant, one may alternatively also pick the gauge 


cj' = 0 . 


(III.27) 


From (III.22) one then finds that the path integral must contain a ghost determinant 
det 


—= Y det As usual, this can be rewritten as a path inte¬ 

gral over a real anticommuting scalar ghost 

R 




d-l 


(III.28) 


This choice may seem more natural but for our purposes it is less useful. The reason is that if 


we set a' = 0, in the Hessian (III.Il) there remains a kinetic term for oj which depends explicitly 


on V, whereas if we set a; = 0 all kinetic operators are independent of V. Since our purpose is 
precisely to avoid singularities due to the appearance of V in the kinetic operators, it is clear 
that for us here the second choice is preferable. 

C. Digression on Einstein-Hilbert gravity 


Since physical gauges are not very familiar, in this section we make a little digression to test 
our procedure in a setting that is better understood. We consider the special case when F and 
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V are constant, in which case (II.3) can be used. In this section we also drop the scalar field 
entirely. The beta functions of the dimensionless couplings A = A/k^ and G = have the 

general form 


dtG = {d- 2)G + BG^ , (III.29) 

dtA = -2A + ^AG + BGA , (III.30) 


The coefficients A and B can be written as A = Ai — r]A 2 , B = Bi — r]B 2 , where rj is the 
anomalous dimension of h^u- 

We begin by considering only the case A = 0 and we focus on the coefficient Bi. We work 
in arbitrary dimension but we are especially interested in the case d = 2 + e. Let us first recall 
the situation with a standard gauge hxing term 

Sgf = ^J d^x^r'^x^.Xu ; , (III.31) 


depending on two parameters a and fd. (The case a = l,/3 = d/2 — 1 has been discussed recently 
in [30].) In the linear parametrization one finds, in the limit d —)• 2, independently of a, 


2(19-38/3 + 13/32) 
3(1-/3)2 


(III.32) 


which for /3 —?• 0 gives B = —38/3, a number that has been found many times in the literature 
|35H40j . The same calculation with the exponential parametrization leads, again independently 
of a, to 


2 (25 - 38/3 + 19/32) 
3(1-/3)2 


which for /3 —)> 0 reproduces the well-known result of 2-dimensional quantum gravity Bi 


(III.33) 
= -50/3 




What is one to make of this discrepancy? We have computed the beta function of G in 
generic dimension d, where it is not universal, and then have taken the limit d —)• 2. In this 
limit G becomes dimensionless and it is generally the case that the one-loop beta functions 
of dimensionless couplings are universal. In fact the limits d —>• 2 do exhibit some degree of 
universality, insofar as they can be shown to be independent of the choice of the cutoff function. 
The choice of parameterization does affect the limit, however. The fact that the limit of the beta 
function does not generally agree with the 2-dimensional result is probably not too surprising, 
since our calculation takes into account all the degrees of fredom of the metric, including the 
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transverse traceless fluctuations, whereas in two dimensions those degrees of freedom do not 
exist, even at a kinematical level. From this point of view it is perhaps surprising that one 
can reproduce the 2-dimensional result at all. We observe that within the approach described 
here (one loop background-field calculation in so called-single metric approximation) the way to 
do so is to use the exponential parametrization and the gauge /3 = 0, whose special feature is 
that it does not involve the conformal degree of freedom. (For the limit it makes no difference 
whether one chooses /3 = 0 for all d, so that the conformal degree of freedom is not involved in 
any dimension, or the de-Donder condition /3 = d/2 — 1, where this is only true in d = 2.) 

Now let us see what happens in the physical gauges. If we use a linear parametrization of 
the metric then the gauge choice = 0, cr' = 0 gives Bi = —38/3. In fact it can be seen that 
it gives the same result for Bi as the standard gauge /3 = 0, a = 0, in any dimension and for 
any value of A. This is because a = 0, (3 = 0 means that we strongly set to zero the quantity 
On the other hand the gauge choice = 0, w = 0 gives Bi = —26/3. We note that this 


coincides with the limit of (III.32) for /3 —>■ oo. 


If we use the exponential parametrization of the metric, the gauge = 0, a; = 0 gives 
Bi = —38/3 while the gauge = 0, a' = 0 gives Bi = —50/3. ^ This calculation confirms the 
previous conclusion that in order to reproduce the 2-dimensional result one has to use a gauge 
condition that does not involve the conformal factor. Since all these results depend critically on 
the form of the ghost terms, we take their consistency as a confirmation of the correctness of 
the procedure. 

We close this section by giving the complete beta functions of A and G in the exponential 
parametrization and in the physical gauges, calculated with the optimized cutoff. In the gauge 
= 0, w = 0 the coefficients of the beta functions are (see [2] for other cutoff types) 

167r(d — 3) 


= 


^2 = — 
Bi = - 

B2 = 


(47r)rf/2r[d/2] 

167r(d — 1) 

(ATT^Gid + 2)T[d/2] 

I67r(d5 _ 4d^ - 9^3 - 48d‘^ QQd + 24 

(47r)'=*/2l2d2((i- l)r[d/ 2 ] 

167r(d3 15^3 58^2 43 


(111.34) 

(111.35) 

(111.36) 

(111.37) 


(47r)«'/2l2d2((i - l)(d + 2)T[d/2] 

We will not attempt to calculate the anomalous dimension here. If we restrict ourselves to the 


® One has to be careful with the order of the limits. Here we are always taking first A —>■ 0 and then d —>■ 2. 
Taking the limits in the opposite order would give, in this case, Bi = —52/3. 
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one-loop terms Ai and Bi there is a Gaussian fixed point with critical exponents 2 and 2 — d 
and a nontrivial fixed point with critical exponents d and d — 2. In d = 3 the non-Gaussian 
fixed point is at A = 0, G = 0.307. In d = 4 it is at A = 0.133, G = 3.35. Two features have 
to be stressed: the reality of the critical exponents and the absence of singularities in the beta 
functions. Both are due to the fact that the beta functions are independent of A. The presence 
of singularities is usually attributed to deficiencies of the truncation and there have been claims 
that the singularity should be replaced by an infrared hxed point [Ml E] . It is interesting to 
find at least one quantization procedure where this troublesome feature is absent even in the 
simple “single-metric” Einstein-Hilbert truncation. ^ 

In the gauge = 0, cr' = 0 the beta functions depend on A (via the kinetic operator of w). 
They are much more complicated and we will only give here their form in d = 4: ® 


Ai 

A 2 




B 2 


3-4A 
7r(3 - 8A) 

-9 20A 

67r(3 - 8A) 

-477 -h 2448A - 33921^ 
247r(3 - 8A)2 
171 - 888A -k 1280A2 
277r(3 - 8A)2 


(111.38) 

(111.39) 

(111.40) 

(111.41) 


We note that in this case there is a singularity at A = 3/8. The non-Gaussian hxed point occurs 
at A = 0.136, G = 2.66 and has complex critical exponents 2.44 ± 0.593z. 


D. Flow equations 


The only effect of the gauge choice is to remove the four degrees of freedom u and ^ from the 
theory and to add a real transverse vector ghost and a real scalar ghost. The case of unimodular 
gravity is essentially identical, except that the scalar ghost is not present. (For recent related 
discussions of unimodular gravity see (Ml US]-) The gauge hxed hessian is 


dxy/g 




fiiA 


-V^ + 


2R 


d(d- 1) 


h 


(^d — l)(^ — 2) 
" 4d2 


F(0)a'(-V2 


cr 


In [2] it was shown that the singularity is avoided also with a “type III” or “spectrally adjusted” cutoff. The 
price one had to pay there was a redefinition of what one means by high and low momentum along the flow. 
Here the same is achieved with a “type I” cutoff, where the distinction between high and low momentum modes 
is conceptually cleaner, being based on the spectrum of a fixed second order differential operator. 

® In d = 3 the non-Gaussian fixed point occurs at A = 0, G = 0.257 and has critical exponents 2.516 and 1. 
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-F'(,^)^W(-V2)(-V2- 




III.42) 


The field a' is invariant under volume-preserving difFeomorphisms, so all three fields are physical. 
As a final simplification we note that defining 


// / I 

a = a + 


2d F'{(j)) 

F(0) 

the gauge fixed hessian becomes diagonal: 


-v^-A 


-V2 


64> 


(III.43) 


dxy/g 




flU 


-V^ + 


2R 


h 


j-'j'fii' ((i — 


1 


+ 4^^ ( -V2 + V"{^) - F"{^)R + 2 


d{d-l) 

d-lF'{^)^ 


Ad? 


F{(j))a"i-V^)a'‘ 


-V^ - 


R 


6(j) 


(III.44) 


d-2 F{(j)) V d - 1 

From here on we proceed as in [26] and for notational simplicity we shall remove the bars from 
the background fields. We choose the cutoff in such a way that the modified inverse propagator 
is identical to (III.44) except for the replacement of —by Pfc(—V^) = —+ Pfc(—V^). We 

I, would amount to a slightly 


note that applying this procedure directly to Eq. III.42 


as m 


different definition of the cutoff. Both procedures seem legitimate, and our choice is dictated 
purely by later convenience. ® 


The flow equations for V and F can be extracted from (1.2) using the formula 

1 


TrIT(-V2) = 


(47r)'^/2 


Q,{W)Bo{-V^) + Qd_,(IT)P2(-V") + ... 


(III.45) 


where the coefficients Qn{W) are given, for n > 0, by Qn{W) = dz z^~^W(z). Since we 

keep at most terms linear in R, we need only the first two terms of the expansion of the heat 
kernel of —V^. We refer to the appendices in |2| for a pedagogical discussion. With standard 
procedure, neglecting in our LPA truncation the anomalous dimension of the scalar held, one 

2-d 

arrives at the flow equations for the dimensionless functions of the dimensionless held (p = k 2 (j)-. 
f{(p) = k‘^~‘^F{(l)) and v{(p) = k~'^V{(f)). 

We shall consider two approximation schemes. As a first case we neglect derivatives of F^ 
with respect to k in the r.h.s. of the flow equation. The analysis of the scaling solutions of the 
resulting equations, and their eigenperturbations, for d = 3 and d = 4 is given in Sections IV 
and V. With the insight obtained in this way, in section VI we shall consider the full equation 


® Neglecting F on the r.h.s. of the flow equation, the fixed point eqnations derived from the two procednres turn 
ont to be the same. 
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where the terms proportional to F are not neglected. This means that we replace dtf —)• 0 in 
the r.h.s. of the fixed point equation (but not of the flow equation, hence also not in the analysis 
of eigenperturbations). The discussion here will be short. In the rest of the paper we consider 
only the cases d = 3 and d = A. In appendix C we shall give the form of the flow equations for 
general dimension. 

If we neglect F, in the r.h.s. the flow equations in d = 3 read 


i 1 , 37 (/ + 4/2)(2 _ 4/// + 5^//) + 

/ ^ ^ ^ ^ 72 - k ‘^ ^ ^ 247r2(/(l+ u")+4/'2)2 


whereas in d = 4 

i) = —Av + ipv' + 


+ 


/ + 3/ 


/2 


167r2 327r2 (3/'2 + /(I + v")) 


r V (/ + 3/'2)(l-3/" + 3i/') + 2/i,« 

/ = -v + ff +^ + f- 


(111.46) 

(111.47) 

(111.48) 

(111.49) 


967r2(3/'2 + /(I + u"))2 
It is interesting to note that if we assume /' = 0, and discard the constant terms in the 


r.h.s, the equation for v reduces to the flow equation (II.1) for the potential in pure scalar 


theory, derived with the same cutoff. In general, shifting the r.h.s. of (II.l), (III.46) or (III.48) 


by a constant results only in a constant shift of the potential of the solution, and is therefore 


immaterial in flat space. The constants in the r.h.s. of (III.47) and (III.49) cannot be discarded 


in the same way because these equations contains / also in the nonlinear part. 


IV. SCALING SOLUTIONS IN d = 3 


A. Analytic solutions 


Let us look for solutions of the system ( III.46|III.47 ). Assuming that v and / are constant, 
there is a unique solution 


1 


n* = 


0.005629 ; 


f* = 


43 


0.06051 


(IV.l) 


187r2 727^2 

In [26] this was called the “Gaussian Matter Fixed Point”, since at this fixed point the EAA 


becomes quadratic in (p. ^ Here we shall call it “FPl”. As in (|II.3|) we change variables to 


V = 


2A 

IQttG 


f = 


1 


IbvrG 


(IV.2) 


^ As pointed out in [46], graviton loops will generally induce derivative interactions in the matter sector, so the 
existence of this Hxed point is limited to a truncation where such interactions are absent. 
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whose beta functions are 


(IV.3) 

(IV.4) 




Ovr 


86 


= -2A - —GA + —G . 


Ovr 


Stt 


In these variables one sees the Gaussian fixed point at (5 = A = 0 (which cannot be seen in the 
other variables) and FPl at 


~ 2 ~ Ott 

A* = — pa 0.04651 ; = — ^ 0.3288 

43 86 


(IV.5) 


The difference with the values given in section III.C is due to the presence of the scalar field 
in the loops. These values are somewhat larger than the ones that had been found for pure 
gravity in three dimensions either in the Einstein-Hilbert truncation, or with a Chern-Simons 
term |22l |23] , or with higher derivative terms [2Tj . 


For this fixed point it is possible to perform an analytical study of the eigenperturbations 
in the infinite dimensional functional space spanned by v and /. The general linear equations 
for the eigenperturbations are of the form (O — X)5w = 0, where SuF = {5v, 5f) and O is the 
corresponding differential operator, which is constructed by substituting into the fixed point 
equations 

v{if) = v^{^) + e5v{ip)e^^ , f{ip) = /*(¥?) + e 6 f{ip)e^^ , (IV.6) 

and expanding to first order e. Linearization around FPl leads to the equations 

r // 

-(A + 3)6v + '^ 6 v' - = 0 

2 bvr^ 

-(A + W+|i/'-g + |4=0 (IV,7) 

which can be studied both analytically and numerically. Imposing that the eigenfunctions 
are even functions of ip leads to a quantization condition for the eigenvalues: A = —9 = 
—3, —2, —1,0,... (with 9 = —A the critical exponent). The relevant perturbations (those with 
0 > 0) are 

6»i = 3, w\ = { 6 v, 6 f)i = (1,0) 

02 = 2, wi = {6v,6f)2= + 

6»3 = 1 , wl = { 6 v, 6 f )3 = ( 0 , 1 ), 
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wi = {Sv,6f)4 = - — 




04 = 0 , 


U >5 = {6v, Sf)5 = ( 0, + (f- 


1 


wl = {6v, 5f)e = (-g^ + (IV. 8 ) 

Note that the eigenspaces of the eigenvalue —1 and 0 are two-dimensional. The asymptotic 
behavior of wq as well as of all the eigenperturbations with eigenvalues A > 0 is power-like. 
We have computed also some eigenperturbations associated to the irrelevant directions but we 
do not find it necessary to give their expressions here. Moreover it can be shown that all 
the eigenperturbations found are not redundant. A check with a numerical eigenperturbation 
analysis has been done using a spectral method based on Chebyshev polynomials and leads to 
the same result. We shall use the numerical approach later on when the analytical approach is 
not possible. 

There is then another analytic solution with constant v but nonconstant /: 

1 __ „ 37 1 


u* =-^ ss 0.005629 ; 


f* = ^ ~ 0.0520678 + 0.25cp‘^ , (IV.9) 

727r^ 4 


which corresponds to 

~ 97r ~ 2 1 

G, = — ^ 0.3821 ; A, = - ps 0.05405 ; ^ = - 

(where we follow the standard terminology of calling ^ the coefficient of the nonminimal coupling 
We shall call this hxed point FP2. We emphasize that although this solution only has 
three nonvanishing couplings, it is a solution of the full functional equations, not of a truncated 
subset. 


The linearized flow equations around FP2 are: 


— (A -|- 3)5v + — 

-(A + i)v+|r- 


(18.V+37) ^ ^ 

(907r^(/9^ -|- 37) 

+ 37) _ , _ 

247r2 (907r2<^2 + 37 ) ) 


(IV.IO) 


We proceed here with numerical methods which we have tested with the previous case. The 
numerical approximation for the eigenfunctions is good for small values of (p but this is enough 
for our purposes. We find four relevant direction. The critical exponents and corresponding 
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FIG. 1: Top: the eigenperturbations 5 v 2 and 6/2 for 9 = 1.815. Bottom: the eigenperturbations Sv 4 and 
6/4 for 9 = 0.516. 

approximated eigenpertubations (rewritten as simple polynomials in a small range aronnd the 
origin) 

64 = 3, w{ = { 6 v, <5/)i = (1,0) 

02 = 1-815, w\ = {5v,5f)2 

03 = 1, wl = {6v,5f)‘i = {Q,l) 

04 = 0.516, w\ = { 5 v, 6 f )4 (IV.ll) 

We show in Fig. [^in different plots the components of W 2 and W 4 . 

Finally there is a fixed point that we shall call FP3: 

1 43 

^;* = ~ 0.005629 ; f, =^-0.0757p^\ (IV.12) 

lovr^ 568 

where Newton’s coupling, defined in terms of /[O], is formally infinite and / is otherwise always 
negative. We shall not give the results for the fluctuations around this fixed point here. 
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B. Search for gravitationally dressed Wilson-Fisher FP 


Finally we look for a solution with nontrivial v and /. We begin by considering finite 
polynomial truncations 


riv "■/ 

v{ip) = Y, ; f{^) = Y • (IV.13) 

2=0 2=0 

We have considered polynomials in of order up to 9 for v and 8 for /. At each order one 
finds a fixed point whose features are close to those of a fixed point in the lower truncation, so 
we assume that these are the trace of a genuine fixed point. For the highest truncation the fixed 
point solution is 


vipp) = 0.00679- 0.0856v?2 + 0.568(^^ + 1.353(/7® + 2.903vp® + 3.390¥5^° 

-7.737(^^2 - 55.00^9^^ - 122.8^9^® + 30.45v9^® (IV.14) 

f{ip) = 0.0625 - 0.0578(/92 - 0.0608^9“^ + 0.00270^9® + 0.318/ + 0.675v9^® 

-2.411/2 _ 20.53<y9i^ - 52.61^9^® . (IV.15) 


This already shows the striking difference between the equations obtained here and the ones of 
pn] . that only admitted potentials unbounded from below. 

The fixed potential has u(0) = 0.00679 and a nontrivial minimum at ~ 0.245 equal to 
V ~ 0.040. Although the fixed point is encouragingly stable as one increases the order of the 
polynomials, the radius of convergence of this series expansion does not extend much further 
than the nontrivial minimum. One can get a slightly improved series by expanding around the 
minimum instead of zero. In this case we obtain 

vipp) = 0.0040+ 0.8877(/-/2 + 2.177(/-/^ + 3.451(/-k)'^ 

-1.653(/ - Kf - 20.03(/ - /® + 12.10(/ - k)^ + 283.9(/ - nf (IV.16) 

-53, 68(/ - Kf - 5592(/ - - 6163(/ - + 12691(/ - . 

f{pp) = 0.0588 - 0.06488(/ - / - 0.05302(/ - + 0.08890(/ - nf + 0.3073(/ - nf 

-1.073(/ - Kf - 5.403(/ - k)® + 17.25(/ - k)^ + 126.0(/ - k)® (IV.17) 

-318.6(/ - nf - 4270(/ - + 1088(/ - , 

where k = = 0.0600, in agreement with the expansion around the origin. This expansion 

reaches slightly larger values of the field than the expansion around zero. 
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It is instructive to compare these approximate solutions to the solution of the equation (II.l) 


which is an approximation of the Wilson-Fisher fixed point in the Local Potential Approximation. 
When compared, the two curves are practically indistinguishable, up to a little beyond the 


minimum of the potential. This can be understood by noting that if /' = 0 the equation (III.46) 


reduces to (II.l). Thus the difference between the two solutions is entirely due to the fact that 
/ is not constant. 

One can calculate the critical exponents in this polynomial approximations by taking the 
eigenvalues of the (n„ + nj) x (n„ + ?^/) matrix of derivatives of the beta functions with respect 
to the couplings. The best estimates give the critical exponents: 3, 1.562, 0.997, —0.588 ± 
0.108i... The first eigenvector points precisely in the direction of Aq, the others are admixtures 
of all couplings. 

These finite truncations can be complemented by numerical studies along the lines of m- 
One tries to solve numerically the fixed point equations u = 0, / = 0 with the initial conditions 
f^(0) = 0 and f'{0) = 0. One can then plot how far the numerical routines can go, as a function 
of the initial conditions u(0) and /(O). Fixed points then typically appear as spikes in this 
graph. We have charted an area around the origin in the u(0)-/(0) plane. There is a large 
region for u(0) < 0 and /(O) > 0 where the solution easily extends up to arbitrarily large (p. 
These solutions all have potentials that are unbounded from below (they behave asymptotically 


as in equation (IV.20) below). For /(O) < 0 there are areas where the behavior look quite 
chaotic. We cannot say much about the system for such initial conditions. The most interesting 
area is a mountainous triangle in the quadrant u(0) > 0 and /(O) > 0, enlarged in Fig. It 
is relatively smooth but the ridges become quite sharp near its vertices and there are distinct 
peaks at the end of each ridge. Two of these can be seen clearly if we cut the graph along the 
line u(0) = l/(187r^), which is common to the two fixed points ( IV.l|V.10 ). Then the two fixed 
points appear as very clear spikes at initial conditions that agree numerically with the values of 
Ao and of the polynomial solutions, as well as the “exact” values. ® The nontrivial fixed point 
can be seen by cutting along the line u(0) = 0.0068, as seen in Fig. One should be careful 
in interpreting such plots. They probably reveal as much about the workings of the numerical 
integration algorithm than about the equations themselves. Nevertheless, the coincidence of the 


* We note in passing that the fixed point FP3, which also lies on the same line, cannot be seen by this technique 
because it corresponds to the initial condition /(O) = 0 where the equations in normal form have a singularity. 
For the same reason this solution also does not show up among the solutions of the polynomial expansion 
around = 0. 
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FIG. 2: Plot of the maximum value of ip reached by the numerical integrator in the most interesting 
region of the initial conditions 0.0055 < ^(O) < 0.0070 and 0.050 < /(O) < 0.065. The jagged appearance 
of the top ridge is a numerical artifact. 




FIG. 3: Left panel: Spike plot as function of /(O) for ri(0) = I/IStt^ « 0.0056. The peak on the 
right is FPl (IV.I), the one on the left is FP2 ( |V.10 ). Right panel: Spike plot as function of /(O) for 
i;(0) = 0.0068. There is no clear spike for intermediate values of r’(O). 


results from polynomial truncations and numerical integration is good evidence for the existence 
of fixed points. Two of these can also be seen analytically, so the main issue is whether the third 
one can be continued to large ip. 

In order to understand the behavior of the solution for large field, let us write the fixed point 
equations v = 0 and / = 0 in normal form: 


/" 


(/ + 4/^) (1 - 37r^(6u - 

37r2/(6u — ifv') 

+ (367rV/' - 727r2/ + 46) (/ + Af'^)) (IV.19) 
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If we assume that for large ip, 1, /, iff and are all negligible with respect to Gv — ipv', then 
the solution would behave asymptotically like 


vr^-^ip^ + Clip + C 2 + 2C3(3^ - ~ CsiSip - 


(IV.20) 


They are analogous to the solutions of the pure scalar flow equation u = 0 with v given by (II.l), 
with initial condition u(0) < 0, or equivalently u"(0) < —1. Such solutions exist for all ip but 
are unbounded from below and are therefore unphysical. 

Another possible asymptotic behavior, which is expected for dimensional reasons, is v{ip) = 
Aip^ + ... and / = Bip'^ + ..., for some constants A and B. This is the behavior that leads 
to the Wilson-Fisher fixed point in the pure scalar case. In order to better approximate the 
solution for large field one can use an expansion in 1/ip. The subleading terms can be calculated 
iteratively from the equations (III.46|lII.47). One finds ® 


, , ,6 1 165 + 1 

v{ip)=Ap> + 
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69127r'^53+8647r^52+277r‘^5 - 3610A 


9'7r2 DGG'K'^Aip^ 24307r‘^A(^® 


^IV..21) 


/(+) = Bip'^ + 
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165 + 1 


+ 
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10206007r6A25(^8 

103687r^53+2167r^52-277r^5+^2^^22^ 


367r2 21GGTT‘^Aip^ GAmir^Aip^ 29 1 60007r6A25y?8 

There are corrections to this behavior which can be obtained by studying the linearized 
differential equations around such an asymptotic solution. One can see that there are four 
linear independent solutions to the linearized equations, but two of them just renormalize the 
constants A and 5 of the previous expansion, while the other two are truly new corrections with 
an essential singularity in the variable 1/ip. We hnd that both 5v and 5f are proportional to 


exp 


<^10 236V® 

1536? 


cosh (-^ip'^ 
\ UV2c 


(IV.23) 


where 6 = 


and c = and differ starting from the power-like term, i.e. in the log 

term inside the exponential. 

Equations! III. 46|III. 47 ) do not have fixed singularities, but they can have movable singular¬ 
ities. At such points ips, both v” and f” become singular; 


v{(p) = {ip - ipsf'^‘^ [^0 + 0{ip - +<()] +Vo + Vl{ip- (fs) + 0{ip- ipsf 
/(+) = (+ - [Bo + 0{ip - + 5 )] + /o + /l(+ - iPs) + 0{p- Psf 


(IV.24) 


In actual calculations we have pushed this expansion to order 1/ip*^. 
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with 


Ao — —8 Bq , Bq — 


V^97r2/o (184 - 97rV^)-12967r4/^-529 


, 36/o-P 

= 6 — fi = —^- 

(fs 18(/?s 


54\/67r3^/^ 

We have not investigated further the equations for other possible singular behaviors. 

If the solution does not end at such a singularity, it should be possible to match the expansion 
for small field to the expansion for large field. In the pure scalar case this method leads to a 
very good overlap of the approximate solutions and to good estimates of critical exponents 
HZllIB]. In the present case this does not happen: the radii of convergence of the expansions 
do not overlap. We have tried to extend the polynomial solutions towards the right by selecting 
suitable Fade approximant. It is then possible to adjust the free parameters A and B of the large 
field expansion to match quite well the small field Fade approximants. We find the best match 
for A ~ 3.41, B « —0.144. Unfortunately this is misleading: neither the Fade approximant not 
the asymptotic expansion are good approximations of the solution near that point. The reason 
is that the Fade approximant of / has a zero at a point (pQ ~ 0.9, and the equation is singular 
for / = 0, as best seen from the normal forms ( IV.18|IV.19 ). 

One can obtain a better approximation of the solution near the zero crossing by using a 
polynomial expansion. The free parameters in this expansion are (pQ, the position of the zero of 


/, and vq = v{po). If one sets f{po) = 0 in (III.46) and (III.47), one obtains v'{po) and f'{po) 
as functions of po and vq. Taking one derivative of (III.46) and (III.47) one determines v"{po) 
and f”{po) as functions of po and vq, and so on. The resulting expansions read 


v(p) = uo+ 0.03377(177.65uo - 1)^^—^ + 0.07810(177.65uo - 1)^^—^^ + .(IV.25) 

‘Po Po 

f(p) = -0.10414^^ +(-0.05548+ 1.300uo + 0.0187v9g)^^-^^^ + ... (IV.26) 

Po Po 

We then integrated numerically the fixed point equations from po towards lower and higher 
values of p. In the direction of decreasing p this can be done all the way to y? = 0 without 
undue difficulties. One then fixes uniquely the initial conditions po and vq by demanding that 
u^(0) = 0 and f'{0) = 0. The resulting numerical solution agrees to high precision with the 
numerical solution determined by the spike plot and also with the polynomial solution. These 
solutions are shown in Fig. The behavior in the region around po is enlarged in FigJ^ The 
initial conditions corresponding to the fixed point are po = 0.92299 and vq = 2.12474. We thus 


in actual calculations we have pushed this to order {p — po)^^- 
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FIG. 4: The solution up to ipQ. Left: v, right: /. Red, continuous curves: the polynomial approximation 
( IV.16|IV.17 ); blue, dashed curves: the Fade approximant; thick black, continuous curves: the numerical 
solution. The potential has a minimum at « 0.245. 


have a reliable numerical solution all the way from y? = 0 to (/? = (/?o- This numerical solution can 
be continued up to y? ~ 35, but it is not at all clear that its behavior reflects the actual solution. 
In the region 1 < (^ < 35, u" reaches a minimum very close to —1 and f passes through a 


zero. Similar clearly unphysical behavior is also observed in the numerical solution of (II.l), so 
we tend to discount it as due to a numerical instability. From scalar theory we also know that 
it is impossible to numerically integrate the fixed point equation from small to large tp. One 
possibility would be to match the asymptotic large-:^ behavior to the expansion around ipQ. If 
one tries to do that, one gets a reasonably good match for v but the match for f gets no better 
than a few percent. This is because the radius of convergence of the expansion around p>Q is 
zero, as one can see by analyzing the coefficients of the expansion, and the asymptotic expansion 
is not good near pQ. The better procedure would be to start at some large p with initial values 
dictated by the asymptotic behavior ( IV.2I||lV.22 ) and to numerically integrate the fixed point 
equation from large to small p. So far our efforts to do this have been unsuccessful. 


V. SCALING SOLUTIONS IN d = 4 


In this section we look for solutions of the system ( III.48|[ni.49 ). First we find that also in 
four dimensions there is a fixed point with constant v and /, which we call again FPl, at 

3 _ „ 41 


u* = 


1287r2 


0.00237 


h = 


7687r2 


0.00541 . 


(V.l) 
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FIG. 5: Enlargement of the solution near Top: v and /. Bottom: their second derivatives. Red 
continuous curve: numerical solution; black dashed curve: the polynomial approximation (IV.25|IV.26|. 


or equivalently at 


A* 


9 

41 


0.219 ; 


G, 


dSvr 

IT 


3.678 . 


(V.2) 


This scaling solution looks like a fixed point of an Einstein-Hilbert theory. If analysed in such 
a two dimensional subspace, then the critical exponents are real and simply given by = 4 
and 02 = 2 with the corresponding eigenvectors in terms of (A, G) components wi = (1,0) and 
W2 = 

In terms of A and G there is also the Gaussian fixed point at (0, 0) with critical exponents 
01 = 2 and 02 = —2 and the same eigenvectors wi = (1,0) and W 2 = (y|^, 1). Actually the 
equations are so simple that one can solve for the RG trajectories even far from the fixed points: 


G{t) 

m 


_ G{to)G, _ 

G{to) + (a - G{to)) e20o-t) 

^ G{to) + (^A(to) - G(to)) 
G{to) + (g, - G{to)) e20o-i) 


(V.3) 


Gomplete RG trajectories (with —oo < t < oo) exist for initial conditions such that G{to) < G*, 
otherwise a Landau pole is present. We note that the separatrix which connects the Gaussian 
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fixed point to FPl is a straight segment, along the direction of W 2 - One can see that the 
corresponding dimensionful couplings G ( 5 e 2 (to t) A ~ Ae for G{to) < (5* reach a 

constant value in the IR. 

Let us recall that the FPl solves the hxed point equations ( III.48|III.49 ). It is therefore 
important to study the eigenperturbations in such a space. For this particular solution the 
linearization leads to the equations 


— (A + 4)(5u + (p5v — 


6 v" 

327r2 


= 0 


Sf" Sv" 

-(A + 2W + ^i/'-^ + ^ = 0 


(V.4) 


whose solutions can be investigated both analytically and numerically. Imposing that the 
eigenfunctions are even functions of (p leads to a quantization condition for the eigenvalues: 
A = —9 = —4, — 2 , 0,... (with 0 = —A the critical exponent). We find 


Oi = 4, w\ = {5v, 5f )i = (1,0) 

6*2 = 2 , wl = { 6 v, 5 f )2 = ( 0 , 1 ) 

03 = 0 , wl = { 6 v, 6 f) 3 = + 


(V.5) 


The dimension of the UV critical surface is 2 without taking into account the marginal direction, 
whose behavior should be further analyzed. 

Since the relevant directions of this fixed point are related to scalar field-independent eigen¬ 
perturbations, if the bare action in the UV is located in the functional space close to it, the flow 
towards the IR may tend to amplify the constant components of v and /. In such a case quan¬ 
tum effects could lead to a reduced field dependence in the effective action, where all quantum 
fluctuations are integrated out. Note that for effective actions with a reduced ip dependence the 
equations of motions in such a model would have solutions closer to be maximally symmetric 
spaces, e.g. de Sitter geometry. 

There is then a nontrivial analytic scaling solution FP2 with constant v but non-constant /: 

3 _, 37 1 


u* = 


1287r2 
which corresponds to 


0.00237, M^) = 


7687r2 6 


-F PS 0.00488-^0.167(^2 , 


(V. 6 ) 


G* = ^ « 4.08 ; ^ ~ 0-243 ; ^ ^ . 


(V.7) 
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The linearized equations for the eigenperturbations around this fixed point read: 


— (A + A)5v + (p5v' 
-{X + 2)5f + ^ 6 f 


(l287rV^ + 37) 
327r2 (3847rV^ + 37) 
(l287rV^ + 37 ) 
327r2 (3847rV^ + 37) 


6 v" = 0 

{ 5 f - 6 v") = 0 . 


(V.8) 


Also here, as in the case with d = 3, we proceed with numerical methods., which give a good 
approximation for the eigenfunctions for small values of ip. We hnd three relevant direction with 
critical exponents: 6 i = 4, 02 = 2, 63 = 1.809. In particular 


01 = 4 , w\ = {6v,6f)i = ( 1 , 0 ) 

02 = 2, wi = {6v,6f )2 = (0,1) 

03 = 1.809, wl = {6v,6f)3 = {ci^ + ip‘^,cif + ip‘^) (V.9) 

where ci^ ~ —0.0028 ad cij ~ 0.53. 

Starting from an initial condition close to this scaling solution a flow towards the infrared 
will enhance such components of v and / in the effective action and in general will develop also 

a non trivial potential. It is not clear if along the flow towards the IR there is a possibility for 

the potential to undergo a dynamical quantum symmetry breaking. We shall return to these 
points elsewhere. 

There is also the analog of the fixed point FP3 that we found in three dimensions. It occurs 
at 

3 41 

~ 0.002374 ; f, = -—cp^ ^-0.09762^\ (V.IO) 

It has the same unappealing properties as in three dimensions and will not be discussed further 
in this section. 

Finally we would like to investigate the hxed point equations in order to search for other 
non trivial solutions, possibly with a non trivial scalar potential. We have used all the methods 
already described in the three-dimensional case. Finite polynomial approximations do not show 
solutions that persist systematically from one truncation to the next, except for FPl and FP2. 
We have also searched for numerical solutions starting from cp = 0, with the initial conditions 
u'(0) = /^(O) = 0. We summarise our findings in a three dimensional spike plot in Fig[^ The 
two spikes that are visible correspond to FPl and FP2. As noted in footnote 8, FP3 cannot be 
seen in this way. No other spikes have been found which lead to a nontrivial potential that is 
bounded from below. 
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FIG. 6: Plot of the maximum value of Lp reached by the numerical integrator depending of the initial 
conditions 0.0023 < ?;(0) < 0.00242 and 0.004 < /(O) < 0.006 . 


We complement this information with an asymptotic analysis. After determining the possible 
asymptotic behaviors of the solutions of the fixed point equations ( III.48[III.49 ), we try to 
construct them integrating numerically from compatible initial conditions in the asymptotic 
region towards the origin. The first four terms of the asymptotic behavior which starts with the 
classical scaling read 


v{p 

fip 


4 1 1 + 125 135A +167r2(l + 125)2 

647r2 2304^7rV^ 589824^5^^4^ ” 

2 15 1 + 125 135^ + 167r2(l + 65-7252) 

^ ^ 2567r2 “ 46087r2^(^2 + l327104^i^]4V 


(V.ll) 

(V.12) 


We have not been able to find new non trivial numerical solutions starting from large values 
of +, but we do find the solutions already known analytically. In doing this we have used 
asymptotic expansion with up to 20 terms for each function. 

We have also tried to investigate solutions for which there exist a value +o where /(+o) = 0 
by analyzing them in terms of a polynomial expansion, similarly to the d = 3 case. The two 
parameters which parametrize the solutions are again and u(+o)- Trying to fix them by 
numerically evolving towards the origin we find that the solution such that u'(0) = f'{0) = 0 
can be reached only for +o oo, in which case one reproduces FPl for small + values. No 
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nontrivial solution with / changing sign seems to exist. Thus all methods point to the same 
conclusion, namely that there are no global scaling solutions in d + 4 beyond the ones we have 
found in closed form. 


VI. SOME RESULTS KEEPING THE F DEPENDENCE ON THE R.H.S. 

We discuss here the flow equations keeping the terms proportional to F in the r.h.s.. We 
proceed in a way similar to the previous analysis of the last two sections, considering the two 
cases in d = 3 and d = 4 dimensions. We find two analytic scaling solutions, which correspond 
clearly to the solutions FPl and FP3 of the previous sections, but the solution FP2 is not present. 
In appendix C we give the general results for d dimensions for the fixed point equations, their 
two sets of analytic solutions and the flow equations. 


A. d = 3 


The fixed point equation for d = 3 reads 

n = .4 ipf (sr + v" + 1) + 10/^^' 

2^ IStt^ bOvT^ (4/'2/ (1 + U')) 

_. l 29^9/ iff (Sf + v" + 1) + Qfv" 

2^ OOtt^ 3607r^/ 1447r2 (4/'^ +/(I + u")) 


iff + 2f^) + 2^/3 ^ (5^/ _ 

307r2/(4(/')" + /(l + ^^"))' 

There is again a solution with constant v and /; 


(VI.l) 


(VI.2) 


u* = 


0.009006 ; 


f* = 


67 


0.07543 


(VI.3) 


457r2 907r2 

Although slightly shifted, this is clearly the analog of FPl and we shall call it again by the same 
name. The linearized flow equations at FPl read 


0 = -(A + 3)Sv + l^Sv' + I (xsf - - gj 

n /^215, a 215 5f 5v" 


(VI.4) 


Studying these equations we find that the scaling solution FPl now has five relevant directions 
and one marginal, with the following approximate critical exponents and eigenperturbations: 


di = 3, 


w\ = {Sv,6f)i = (1,0) 
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(VI.5) 


02 = 2, 

215’ 
04 = 1, 


05 = ^, 

^ 215’ 


w\ = iSv,6f)2 = - 


8 2 67 \ 

+ ^ ’“4^j 


Wig = { 6 v, 6 f )3 = ( -^,1 


277r2 

36 


i / 362 16 2 , 4 8308 67 2A 

w, = (fo,s/)4=(^j^-^v j 

= (<^^) ^f)5 = 


( 200049 


36 


17741 


68 


Vll9961407r2 377^ ’ 3503047r2 123 


= 0 , 


w\ = { 6 v, 6 f)e = - 


1810 1810 2 40 4 6 

42937r6 4777r4‘^ ~ 9 ^^^ ’ 

74935 41540X 2 335 


25 7587r6 42937r41627r2'^ 

The two pairs of degenerate solutions with critical exponents 1 and 0 that we found in section 
IV have split, with one in each pair becoming slightly more relevant. 

There is no solution of ( VI.1[VI.2 ) that could be identified with FP2. The fixed point FP3, 
however, is still present exactly in the same position: 


1 


n* = 


568^ 


(VI.6) 


187r2 ’ 

The eigenperturbations around FP3 are different. We have analyzed them numerically but we 
do not find it very useful to show them here. 


B. d = 4 


In four dimensions the fixed point equations are 


n = -4 + , , ^ _ (t>f {Qf + v" + 1) + Qfv" 

^ ^ 327r2 1927r2/ 1927r2 (3/'2 + / (1 + n")) 

n = - 2 /-+ f + v" + '^) + ^fv'' 

J + + 11527^2 2887r2/ 3847r2 (3/'2 + /(I + r;")) 

iff + n {2f + + 6//^ - 2^fff") 

647r2/ (3/'2 + /(I + v")f 

The analytic solution FPl is given by 




5 

1287r2 


0.003958 ; 


f* 


157 

23047r2 


0.006904 . 


(VI.7) 


(VI.8) 


(VI.9) 


Linearizing the flow equation given in 


Eqs. ( |C.8|C.9D for d 


4 around FPl we find: 


0 = 
0 = 


-(A + 4:)6v + (p 6 v' + ^ {X5f - (p 6 f) 


Sv" 
~ 3^ 
6 v" 
96'7r2 


(VI.IO) 
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Studying analytically these equations we find that FPl has three relevant and one marginal 
direction: 


01 = 4, 

314 

02 = — ~ 2.553, 


03 = 2, 

68 

04 = ~ 0.553, 


05 = 0, 


w{ = {Sv,Sf)i 
w\ = 5f)2 

w\ = {5v, J /)3 

wl = {6v, 6f)4 

wl = {5v, 6f)4 


( 1 , 0 ) 



/_ 29 2 157 \ 

\ 5447r^ ^ ’ 32647r^ ) 

/ 56913 72 2 17741 ^ 

V30943527r2 “ ’ “3503047r2 ^ 

( 87 87 _2 ,A 215 
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+ 


(VI.ll) 
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The fixed point FP2 is absent but FP3 is again present in the same position: u* = 
f*{f) = 11 1^^® relevant directions, with critical exponents 4, 2.104, 1.574, 0.2475. 

The eigenvector corresponding to eigenvalue 4 has components {5v, 5f) = (1,0), the others have 
been studied numerically. We find that if 6 v is chosen to be increasing for large (j), as is required 
by physical considerations, then 5f is negative and decreasing for large field. Given that / is 
already negative, this is again physically undesirable. 

Whether the true scaling solutions of the theory (assuming that any exists) resemble more 
this case or the one described in section V, is something that will have to be studied taking into 
account the constraints inposed by the split-symmetry Ward identities. 


VII. CONCLUSIONS 

In this work we have looked at ways to improve the flow equations for the functions / and v 
characterizing the EAA for scalar-tensor gravity, truncated at two derivatives. We have argued 
in favor of using the exponential parametrization of the metric fluctuation, which is more natural 
if one takes into account the fact that the metric carries a nonlinear realization of the linear 
group, together with a physical gauge fixing that simply removes the gauge degrees of freedom 
from the calculation. 

As usual the EAA has been truncated: aside from the gauge fixing, ghost and cutoff terms, 
we have assumed that the EAA is a function of the full metric and of the full scalar field (/>, 
and not of and 4' and 64 separately. This is known as a “single field approximation” 
and its limitations have come under close scrutiny recently [12( 1481150) . Within the single metric 
approximation but in contrast to much of the previous literature, here we have retained infinitely 
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many couplings, namely those pertaining to the first order of the derivative expansion in the 
metric, and the local potential approximation in the scalar. This expansion is well motivated 
in applications to critical phenomena, less so when studying the UV limit, but this approach is 
still more systematic than additionally truncating on powers of the scalar field. The consistency 
of the truncation should be tested by enlarging it: in the scalar sector the next step would be 
to take into account the wave function renormalization of (p, while in the gravity sector it would 
consist of including all four-derivative terms, or perhaps a truncation of the form F{R,4>), as 
discussed earlier in m- 

We have considered two approximation schemes, one where we dropped terms proportional 
to F in the r.h.s. of the equations and a second one where we keep them. This is similar to, but 
obviously more general than, neglecting G on the r.h.s. of the flow equations for the Einstein- 
Hilbert truncation. Let us stress that the anomalous dimension of the quantum fluctuations of 
the metric is not directly related to /. An approximate computation of the anomalous dimension 
can be obtained from the flow of the two point function of the fluctuation held, along the lines 


The new equations do not have a singularity for u = /, as previous equations did. One 
consequence is that in the two dimensional Einstein-Hilbert subspace (where v and / are con¬ 
stant) there is no singularity for A = 1/2, nor for any other value of A. We expect that this 
feature will be useful to study the IR limit of the theory. Another consequence is that in three 
dimensions the new equations admit, at least neglecting F and for sufficiently small helds, a 
nontrivial hxed point that resembles very closely the Wilson-Fisher hxed point of scalar theory 
in the LPA approximation. We have traced this similarity to the fact that the flow equation for 


V, eq.(III.46), reduces, for constant /, to the flow equation of scalar theory (II.1). Recall that, 


in contrast, the “old” equations admitted only solutions with negative Taylor coefficients 
The function / starts positive at (/? = 0 and decreases monotonically. The solution exists at least 
up to the point (po where / crosses zero. At this point the new equations have a singularity, but 
we have shown by analytic methods that the singularity can be traversed. We have not been 
able so far to join the expansion around ipo to the expansion around infinity. This may be a 
technical issue that could be overcome by using more powerful numerical techniques, but there 
may also be deeper issues due to the vanishing of many terms in the Hessian when F = 0. 

In addition to this tentative gravitationally dressed Wilson-Fisher hxed point in three di¬ 
mensions, we have also shown that, in any dimension, the equations where one neglects the 
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dependence of F in the r.h.s. admit three other simple fixed points: a “Gaussian Matter Fixed 
Point” FPl with constant v and /, a fixed point FP2 with v constant and / = /o + with /o 
and /i positive (in the interesting range of dimensions) and a fixed point FP3 with v constant 
and / = with /i negative. We emphasize that although these fixed points have only a few 
nonvanishing couplings, they are not just the solution of a polynomial truncation but of the full 
fixed point equations for the functions v and /. In three dimensions FP2 is characterised by four 
relevant eigenperturbations while FPl has four relevant and two marginal eigenperturbations. 
In four dimensions we have shown that FPl and FP2 have only three relevant (or marginal) 
deformations, making them good candidates for asymptotically safe models. We also note that 
the critical exponents we find are not related to redundant operators, in contrast to what was 
seen for the flow equation of f{R) gravity in |in) . In any dimension, the fixed points FPl and 
FP3 are also present, with little or no change in their properties, when the terms proportional 
to F are retained in the r.h.s.. Unfortunately the physically interesting fixed point FP2 disap¬ 
pears in this case. We do not know at present whether this reflects a genuine property of the 
underlying physics. 

One may be worried by this and more generally by the different properties of the equations 
that we have found here, compared e.g. to those of |26j or to the flow equation in the physical 
gauge = < 7 ' = 0 (as discussed in the end of section III.B). The situation is similar to that of 
f{R) gravity, where different equations turned out to have rather different solution spaces [8]. 
The answer seems to be that some approximations are too drastic: all equations are good enough 
to find the fixed point within finite dimensional truncations, but the study of its properties in 
an infinite dimensional function space is more delicate and requires better approximations. It 
has been shown in |12) that pathological features resembling those encountered in f{R) gravity 
can be artificially induced even in pure scalar theory by an improper use of the background 
field method. In particular, one should pay close attention to the violation of split symmetry 
, which, at linear level, amounts to dcj) ^ 5(j) + d'ljj, (p ^ (p — 5ip in the scalar sector and 
h^i/ —>■ h^u + a; —)> w + du, g^u —>■ 9 ^lu — in the gravitational sector. While an 

investigation of this point will be necessary, it seems that the equations derived here are already 
powerful enough to discover at least some of the scaling solutions in the theory. This may be 
a hint that, within the single-field approximation, the use of the exponential parameterization 
and of the unimodular physical gauge is to be preferred. 

Let us note that unimodular gravity corresponds to the case where the conformal fluctuations 


(A. 
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h are completely absent. For such a theory, in the single field approximation for the average 
effective action, the flow equation is obtained from the one of full gravity in the unimodular 
gauge by removing the corresponding ghost contribution, which is a constant term in both 
equations for v and /. Since the running of / does not depend on v but only on its derivatives, 
this is consistent with the fact that in unimodular gravity the constant term in the potential is 
an integration constant. In this framework any constant value of v at the fixed point would not 
contain any physical information. 

There are several obvious extensions of the truncation that we plan to return to in the 
future. Also, we have focused here mainly on the mathematical properties of the system of flow 
equations, but ultimately one is interested in physical applications. In this regard we observe 
that the fixed point FP2 in d = 4 has the properties that were discussed in m as prerequisites 
for the construction of interesting cosmological models. With the linearized perturbations given 
here and with numerical integration of the flow equation it will be possible to analyze in detail 
several scenarios. 

Note added. While this paper was being considered for publication, a work by Borchardt 
and Knorr |56] appeared where they describe the use of pseudo-spectral methods to solve linear 
and non-linear ODEs. They describe a numerical solution of the full fixed point equations in 
three dimensions ( VI.l|VlT2 ). In contrast to the solution of the equations ( III.46|III.47 ) that we 
described in section IV.B, this solution has a growing / that never crosses zero. 


Appendix A: Multiplicative quantum-background split 


The metric is generally thought of as a tensor, hence as a “linear” object. This can be 
misleading. In General Relativity the metric is subject to the constraint of being non-degenerate 
and of having a hxed signature. These constraints define a nonlinear subspace in the space of 
symmetric tensors, and the metric is assumed to lie in this subspace. To understand this 
subspace better, observe that a metric can be defined by giving a frame and declaring it to be 
orthonormal. This defines a surjective map from the set of all frames, which is in one-to-one 
correspondence with the group GL{d), to the set of metrics. Two frames that differ by the action 
of a (pseudo)-orthogonal transformation define the same metric. Hence the set of all metrics 
with signature (p, q) at a point can be identified with the coset space GL{d)/0{p, q), which is an 
open subset in the space of symmetric tensors. It is not a priori clear whether the gravitational 
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path integral should be restricted to this subset or could include also degenerate metrics or even 
metrics with other signatures [5l]. We assume here the former point of view. Then, it is not 
clear, when one employs the standard linear background-fluctuation split 

9 ^iu — 9^u T , (A.l) 

how to implement the requirement that the quantum fluctuations should stay within the space 
of metrics of a given signature. 

From this point of view it is more natural to use a multiplicative version of the background 
field method, introduced previously in ([55]). Any metric of given signature in d dimensions 
can be written as 


9iJ-u — f^O u9po 


(A.2) 


where g is a fixed metric of the same signature as g and 9^ n is a field with values in the group 


GL{d). The multiplicative quantum-background split (A.2) has an inherent arbitrariness that 
consists of local GL((i)-transformations 


= A pi{x) ; g'pAx) = p{x)A>^a{x)gaf3ix) ■ 


(A.3) 


This is a multiplicative version of the “split symmetry” g^^ —)■ g^i, + hpy —)■ hp^, — epy of the 


linear split (A.l) and is mathematically similar to a gauge invariance. Note that when gp^, = 5py 
the field 0 can be viewed as an orthonormal frame. (We discuss here the Euclidean case. Aside 
from global issues, other signatures can be treated in the same way). Every gp^ can be regarded 
as a split-symmetry-transform of 5p^, so there is a subgroup of GL{d) that is conjugate to 0{d) 
and leaves g invariant. We shall refer to it as the stabilizer of g. In a neighborhood of the 
identity, one can write 6 = , where X^p is a generic matrix. Let Yp^, = gppXPy. For any g 

and g, it is possible, by means of a split-symmetry transformation belonging to the stabilizer of 
g, to choose the matrix X in such a way that Y is symmetric. Then it is easy to check that 


g = e^"ge^ = g = ge^^ 


(A.4) 


Putting 2X = /i, this is our starting point (III.l) 


There, we were only interested in a perturbative evaluation of the path integral and therefore 9 was subsequently 
expanded additively in small fluctuations around a background 6 . Here we shall not do this. 


35 






Appendix B: Jacobian for the functional measure 


In this paper we have assumed that the functional integral over metrics has a simple measure 
when the field is defined by the exponential parametrization. An alternative point of view 
is that the measure is simple when is defined by the usual linear parametrization. In this 
case there would be a nontrivial Jacobian, which we calculate in this Appendix and, for the 
truncation we have chosen, does not contribute to the flow equations. 

On differentiating the matrix g = ge^ we obtain at each point of space-time 

OO - 

71=0 


where 


n—1 


n—1 


d{h^) = h^{dh)h^-^-^ = Y[h ’'^^ 


-l-k' 


dh 


k=0 


k=0 


(B.2) 


The last equality can be understood easily if we make explicit use of indices (repeated indices 
are implicitly summed): 


iXdhY)ij = XiadhabY^,, = XiaY[,dhab = [X ® Y^] . dKb = ([X ® Y^] dh).. (B.3) 

Therefore we have for the functional integration measure 


[dfif] = Y)et[gM][dh] 


(B.4) 


where 


OO ^ n—1 


= E E ® ('"b 


n—1—fc 


n\ 

n=0 fc=0 


(B.5) 


We can consider a representation of the determinant based on functional integration over 
Grassmann fields: 


Det[M] = J 


(B. 6 ) 


We note that within the Functional RG approach we can employ a coarse-graining also on 
the 9 and 9 fields which can be considered as quantum fluctuations. One can therefore see that 
in the background field approach for a single metric truncation the flow is not affected by the 


extra term in the action given by (B. 6 ). Indeed if /i 7 ^ 0 this is an interaction term which we may 


choose to not include in the lowest order truncation. This measure would affect the computation 
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of the anomalous dimension of the quantum metric. Moreover it would come into play when 
adopting more complicated (bimetric) truncations. 

Alternatively one can compute directly the determinant of M (x) at each point of space-time 
if the matrix h can be diagonalized by a similarity transformation. Unfortunately in general 
h is not symmetric since h = g~^Y where Y is symmetric. Let us note that for the case of 
interest to us, like a maximally symmetric Euclidean background, we can eventually choose a 
coordinate system for which h is symmetric (gauge fixed background). On general grounds if 
we can diagonalise h by a similarity transformation so that h = RA.R~^ with A = diagjAj} a 
diagonal matrix, then M{x) is similar to the diagonal matrix 


n—1 


n=0 k=0 


A'" (g) A' 


n—l—k 


(B.7) 


which has the same determinant of M{x). One finds 


oo ^ n—1 


^ oo - \n _ \n 

detAfw = nE;)E^Lr"‘=nE;^(«^r‘)+nEstvt^ 


n! 

ij n=0 k=0 


i i<j 


\ . \ \ 2 
gA, _ 


= e 


i n=0 
,tr/i 


n=0 


This leads to an explicit form of the Jacobian in the functional integral. 


\, \ . \ 2 


(B.8) 


Appendix C: Some results in d dimensions 


First we consider the scheme where we purposefully neglect the dependence on (the dimen¬ 
sionful) F present inside the trace through the Rk term. lin arbitrary dimension d the flow 
equations are 


i, = -,lv + l(d - 2 )4,,/ + |(d - l)(d - 2) - 


(Cl) 


and 


. , , 1 , , (d^-Ad‘^-9d^-48d^ + 60d + 2A) 

f -{2-d)f + -{d-2)<Pf -c, 2 M{d-l) + 

dcd {2{d-l)f^ + {d-2)f) {2{d - l)r^ + {d - 2)f) {{d - 2)ff"+2f'^) 

I2 2{d-l)f'^ + {d-2)f{v" + l) {2{d-l)f'^+ {d-2)f{v"+ l)f 

where Cd such that = (47r)‘^/^r((i/2 -|- 1). There are three particular scaling solution of the 
fixed point equations. The fixed point FPl for constant v and / which is given by 


V* = Cd- 


{d-l){d-2) 


f* — C-d 


{d^ - - 7d^ - 50^2 + 60d + 24) 

24d(d- l)(d-2) 


(C.3) 
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the fixed point FP2 given by 


{d-l){d-2) ^ ^ , 

-2(i-’ 


where 

/o = -Cd 

and the fixed point FP3, with 


2{d-l) 

{d^ - 4^4 - 7d^ - 50^2 + 84d + 24) 
2Ad{d - l){d - 2) 




(C.4) 


(C.5) 


1 (d-2)(d-l) , (d-2)(d5_4d4_7d3_50d2 + 60d + 24) ^ 

- Td -■ ~ 8(d - 1) (d<' - - 7# - 44^2 + 72d + 2if ^ 


Now we consider the full flow equations. Taking into account the relations between dimen¬ 
sionless and dimensionful quantitites for the terms which are induced by the cutoff dependence 
on F and F', 

d — 2 j^i-d/2pi = f' ~ j' _ d — 2. ^^ 


k^-dF = f + (d-2)f-'^^f. 


-Vf 


(C.7) 


2 " 2 

the flow equations (not the fixed point equations) are more complicated and on the r.h.s. both 


/, /' do appear. Starting again from Eq. (III.44) they are given by 


1 , {d-l)(d^-d-3) {d-2){d+l)hf-{d-2)ipf) 

V = -dv +-{d-2)ipv +Cd- - TV^ - - + Cd - ^ ^ - - (C.8) 


Fcd- 


d + 2 ' 4(d + 2)/ 

2(d2-4)/+(d-2)(l+u") (2/-6/-(d-2)<^/')+4(d-l)/' ( 2 /'+ (d-l)/'-(d-2)(^r 
2{d + 2) {2{d - 1) {ff + {d- 2)f (1 + u'O) 


/ = 


{2-d)f+^{d-2)ct>r-Cd 


d^ - 2d^ - 15d^ - 46^3 + + ged - 24 


-Cd 


l2{d + 2) {d-l)d 

{d^-nd?-md?+M+A3) (2f-{d-2)(pf 


(C.9) 


- Cd {(d - 2)ff + 2/2) X 


48(d- l)d(d + 2)/ 

(d - 2)(d + 2)f + (d - l)/'2 f (d - 2W-2f) + 2(d-l)//' ((2 - d)^r+ (d + 2)f + 2f 


+Cd- 


(d + 2)/ (2(d - 1) iff + (d - 2)f (v" + 1)) 
f (id-2)^ (4(d - 1)/" + (d - 2) iv" + 1)) - 8(d - 1)/') + 2(d - 2) ( d fv" - f f' + 1 


24 (2(d-l)(/0V(d-2)/(u" + l)) 

On setting f = f' = 0 one obtains the fixed point equations. Now one finds only two analytic 
solutions. The first, which we denote again FPl, has constant v and / and is given by 

-d® -f 2d^ + 15d^ -f 46d3 - 38d‘^ - 96d + 24 


vl = Cd (d — 4:+ -h ^ ) , / 

* V d + 2 dj ■'* 


= Cd- 


12(d-2)(d-l)d(d-k2) 


(C.IO) 
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The second is identical to (C.6) and is therefore called again FP3. We note that in FP3 the 


second derivative of / is always negative, for dimensions of interest. The fixed point FP2 has no 
analogue among solutions of this equation, at least among those for which / is at most quadratic 
in (j). It is possible that its role is taken by a more complicated solution that we have not looked 
for. Let us also stress that since the cutoff is dependent on functions of F{(f>) and F'{^) the 
split symmetry Ward identities will modify the flow equations and the fixed points could change 
again. Their role will be the subject of further investigations. 

Finally we consider the equations derived from an Hessian which is not diagonalised, i.e. 


starting from Eq. (III.42), instead from Eq. (III.44) while retaining the complete F dependence. 
Eor this case the flow equation for v is the same while the flow equation for / is obtained from 


(C.ll) 


Eq. (C.9) by adding to the right hand side 

^ _ /' 2(/ /' - /'/) -id-2) iff + ipff" - ipf^) 

^ ''''(d + 2)/ 2(d-l)/'2 + (d-2)/(l + u") 

The fixed points EPl, EP3 given in ( C.10|C.6 ) are fixed points for this flow too, and again there 
is no analog of EP2. This confirms the expectation that the diagonalization process is a rather 
mild change in the coarse-graining scheme. 
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